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Neutrino Masses from Broken R-Parity 
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a Instituto Superior Tecnico, Departamento de Fisica 
A. Rovisco Pais 1, 1049-001 Lisboa, Portugal 

We review models where R-parity is broken, either spontaneously or explicitly. In this last case we consider 
the situation where R-parity is broken via bilinear terms in the superpotential. We show that although at tree 
level only one neutrino gets mass, at one-loop level all three neutrinos became massive. We study the conditions 
under which bimaximal mixing can be achieved and show that the masses can be in the correct ranges needed for 
solving the atmospheric and solar neutrino problems. 
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1. Introduction 

In the past most discussions of supersymmetric 
(SUSY) |l],|| phenomenology assumed R-parity 
(Rp) conservation where, 



R P = (-1) 



2J+3B+L 



(1) 



This implies that SUSY particles are pair pro- 
duced, every SUSY particle decays into another 
SUSY particle and that there is a LSP that it 
is stable. But this is just an ad hoc assump- 
tion without a deep justification. In this talk we 
will review how Rp can be broken, either spon- 
taneously or explicitly, and discuss the most im- 
portant features of these models || . We will also 
describe recent results on one-loop generated 
masses and mixings in the context of a model that 
is a minimal extension of the minimal extension 
of the MSSM-GUT § in which R P Violation 
(RPV) is introduced via a bilinear term in the 
MSSM superpotential ~ 



2. Spontaneously Broken R-Parity 

2.1. The Original Proposal 

In the original proposal Q the content was just 
the MSSM and the breaking was induced by 



\Vt) = vl 



(2) 



The problem with this model was that the Ma- 
joron J coupled to Z° with gauge strength and 
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Table 1 

Lepton number assignments. 
Field L l~ c V c S others 

Lepton #1-1-11 



therefore the decay Z° — > p^J contributed to the 
invisible Z width the equivalent of half a (light) 
neutrino family. After LEP I this was excluded. 

2.2. A Viable Model for SBRP 

The way to avoid the previous difficulty is to 
enlarge the model and make J mostly out of 
isosinglets. This was proposed by Masiero and 
Valle f|. The content is the MSSM plus a few 
Isosinglet Superficlds that carry lepton number, 

v<i ee (1,0,-1) ; Si = (1,0,1) 5 * = (1,0,0) (3) 

The model is defined by the superpotential ||,|l0| , 

W = h u u c QH u + h d d c QH d + h e e c LH d 

+ {h H u H d - /i 2 )$ 

+h v v c LH u + h<$>v c S 

where the lepton number assignments are shown 
in Table [lj The spontaneous breaking of R parity 
and lepton number is driven by |10| 



\VRj 



v s 



{Z>t) 



(4) 



The electroweak breaking and fermion masses 
arise from 



(H u ) = v u 



(H d ) = v d 



(5) 



2 



with v 2 — v\ + v d fixed by the W mass. The 
Majoron is given by the imaginary part of 



Vv> 



■{v u H u - v d H d ) 



V " T - y^r + (6) 



where V 



Since the Majoron is 



mainly an SU(2) <g> U(l) singlet it does not con- 
tribute to the invisible Z° decay width. 

2.3. Some Results on SBRP 

The SBRP model has been extensively stud- 
ied. The implications for accelerator jll| and 
non-accelerator |L2[ physics have been presented 
before and we will not discuss them here ||. As 
in this talk we are concerned with the neutrino 
properties in the context of Rp models we will 
only review here the neutrino results. 

• Neutrinos have mass 

Neutrinos are massless at Lagrangian level 
but get mass from the mixing with neutrali- 
nosQL!]. In the SBRP model it is possible 
to have non zero masses for two neutrinos 



• Neutrinos mix 

The mixing is related to the the coupling 
matrix h Vi . . This matrix has to be non di- 
agonal in generation space to allow 



J 



(7) 



and therefore evading JL4J the Critical Den- 
sity Argument against v's in the MeV range. 

• Avoiding BBN constraints on the m VT 

In the SM BBN arguments |Q rule out v T 
masses in the range 



0.5 MeV < m„ T < 35MeV 



(8) 



We have shown [|16| that SBRP models can 
evade that constraint due to new annihila- 
tion channels 



v r v T — > J J 



(9) 



3. Explicitly Broken R-Parity 

The most general superpotential W with the 
particle content of the MSSM is given by 



W = Wmssm + Wft 
where 

Wmssm = £ a b 



(10) 



and 

W R 



h%QtU 3 H b u + h%Q b D 3 H a d 
Vh%L b R 3 H a d - vHSHt] (11) 



+\ ljk D l D J U k +e ab e l L?H b u (12) 

where i,j — 1,2,3 are generation indices, a, b = 
1,2 are SU{2) indices. 

The set of soft supersymmetry breaking terms are 



V 



Oft _ T/So/t 



MSSM 



soft 



(13) 



m z H HTH a d 

3 



m H u H u* H u 



i=l 



+A%L\R 3 H a d - B^H a d H b u + h.c 



(14) 



and 



T rSO ft 

V n =E ab 



+ A'(. I), !),( ), + e ab BieiLfH* + fc.c(15) 

The bilinear Rp violating term cannot be elim- 
inated by superfield redefinition as sometimes it 
is claimed. To show this we consider the case JlTj 
where all the trilinear couplings in Eq. ( |l2| ) are 
zero and for simplicity we take e± = 62 = 0. Then 
the superpotential is 

W = e ab \h t Q a 3 U 3 H b u + h b Q b 3 D 3 H a d + h T L b 3 R 3 H a d 
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Consider now the rotation 



H' d = 



nH d - e 3 L 3 



VV 2 + ' 



L' 



e 3 H d + fiL 3 

Vv 2 + 4 



(16) 



(17) 



In the new basis 



W = h t Q 3 U 3 H u + h b ^Q 3 D 3 H' d + h T L' 3 R 3 H' d 
A* 

e 3. 



-n'H' d H u + h b ^Q 3 D 3 L' 3 
A 4 



where 



(18) 



(19) 



But the soft terms, 



V, 



soft 



Mi\L 3 



H d \ LL d\ T^L 3 

B [iH d H u - B 2 e 3 L 3 H u + h.c. 



(20) 



in the rotated basis is 



V, 



soft- 



i 2 H u 2 +M? t\ , „ m 2 H e 2 +M? u 2 ~ „ 

' / '-'3 £ | jji 1 2^_ n d i L 3 r | jj 1 2 



,/2 



,/2 



H d H u J2 \ m H d ~ M L 3 ) L 3 H d 

A* 



/'' 



-^E(B 2 -B)L' 3 H u + h.c. 
A* 



(21) 



The last two terms violate Rp and induce a non- 
zero VEV for the r sneutrino field in the rotated 
basis (y' T ) = v' 3 /\/2, where 



[i! 2 m\ 



(v^Am 2 + n'v 2 AB) 



(22) 



and 



Am 2 = m 2 



Hi 

3h 



Mi 



8tt 2 

AB = B 2 - B 



'' 1 m 2 Hl +M 2 +M 2 D +A 2 D )\n 



8^ 2 



b Anln 



M GUT 
mz 



M GUT 

mz 



(23) 



4. Bilinear R-Parity Violation 

4.1. The Model 

The superpotential W for the bilinear Rp vio- 
lation model is given by 



W = e al 



-fiH d H u + eiL^Hl 



(24) 



where i,j = 1,2,3 are generation indices, a,b = 
1,2 are SU(2) indices. The set of soft supersym- 
metry breaking terms are 

v soft = M^ 2 QTQ a 3 + M% 2 uiUj + u'^iri), 

+Mi j2 LrL° + MfR*R 3 + m 2 Hd HTH a d 



Sab 



+A l £L l lR j H d l -BfiH d l Hu + BiCiL^H, 



(25) 

The electroweak symmetry is broken when the 
VEVS of the two Higgs doublets Hd and H u , and 
the sneutrinos. 



H. 



Hd 

^[xi + Vu + ifl] 
72& 



■^mlvF + Vi + w! 



The gauge bosons W and Z acquire masses 



12 2 

j9 v 



4 



(ci 2 + g' 2 )v 2 



(26) 
(27) 

(28) 
(29) 



where 

v 2 = v 2 d + v 2 u + vj + v\ + v 2 3 = (246 GeV) 2 (30) 

We introduce the following notation in spherical 
coordinates: 

v d = v sin 9 1 sin 9 2 sin 3 cos (3 
v u = v sin 9\ sin 6 2 sin 9 3 sin (3 
Vi = v sin 9\ sin 9 2 cos 9 3 
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v% = v sin U\ cos 62 

V3 — V COS 6\ 

which preserves the MSSM definition tan/? = 
v u /vd- The angles 0$ are equal to it/2 in the 
MSSM limit. 

The full scalar potential may be written as 



E 



dW 



dz 



+ V D + V soft + V RC 



(31) 



where Zi denotes any one of the scalar fields in 
the theory, Vn are the usual D-terms, V so f t the 
SUSY soft breaking terms, and Vrc are the one- 
loop radiative corrections. 

In writing Vjic we use the diagrammatic method 
and find the minimization conditions by correct- 
ing to one-loop the tadpole equations. This 
method has advantages with respect to the ef- 
fective potential when we calculate the one-loop 
corrected scalar masses. The scalar potential con- 
tains linear terms 



Vn near — td^d + *« CT u + = *a°a ) 



where we have introduced the notation 



(-0 J ,.R ,,R ,,R\ 



(32) 



(33) 



and a = d, it, 1, 2, 3. The one loop tadpoles are 



t — t' 

L a b r 



St 



MS 



T a {Q) 



where T^{Q) = -St™ 
one-loop tadpoles. 



(34) 

T a (Q) are the finite 



4.2. Main Features 

The e-model is a one (three) parameter(s) gen- 
eralization of the MSSM. It can be thought as an 
effective model showing the more important fea- 
tures of the SBRP-model [[[o| at the weak scale. 
The mass matrices, charged and neutral currents, 
are similar to the SBRP-model if we identify 



vrK 



(35) 



The Rp violating parameters e$ and Vi violate 
lepton number, inducing a non-zero mass for only 
one neutrino, which could be considered to be the 



CO 1 
CO 1 



0.95 




0.9 



|v 3 | (GeV) 

Figure 1. Ratio of the lightest CP-even Higgs 
boson mass rrih in the e-model and in the MSSM 
as a function of V3. 



the v T . The v e and remain massless in first 
approximation. As we will explain below, they 
acquire masses from supersymmetric loops 
that are typically smaller than the tree level mass. 

The model has the MSSM as a limit. This can 
be illustrated in Figure [l] where we show the ratio 
of the lightest CP-even Higgs boson mass mu in 
the e-model and in the MSSM as a function of V3. 
Many other results concerning this model and the 
implications for physics at the accelerators can be 
found in ref. 

5. Radiative Breaking 

5.1. Radiative Breaking in the e model: 
The minimal case 

At Q = Mqut we assume the standard mini- 
mal supergravity unifications assumptions, 

A t = A b = A T = A , 

B = B 2 = A - 1 , 



> 



5 



M% = Ml = Ml = ml , 

M 3 = M 2 = Mi = M 1/a 

In order to determine the values of the Yukawa 
couplings and of the soft breaking scalar masses 
at low energies we first run the RGE's from the 
unification scale Mqut ~ 10 16 GeV down to the 
weak scale. We randomly give values at the uni- 
fication scale for the parameters of the theory. 
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h ^Gur/ 47T 
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A/m 

MG(/t/ to 
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3 
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Mi/ 2 /m 
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io- 2 


< 


e i GUT/ m 
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10 



(36) 



The values of hf UT , hf T , hf UT are defined in 
such a way that we get the charged lepton masses 
correctly. As the charginos mix with the leptons, 
through a mass matrix given by 



M c = 



M c 
D 



A 

M L 



(37) 



where Mc is the usual MSSM chargino mass ma- 
trix, 



Mc = 



M 



x/2 



gvd 



(38) 



Ml is the lepton mass matrix, that we consider 
diagonal, 

l 



Mr = 



^h Ell v d 











-^h E22 v d 











^h Ea3 v d 



■/2 



(39) 



and A and B are matrices that are non zero due 
to the violation of Rp and are given by 



(40) 



We used an iterative procedure to accomplish 
that the three lightest eigenvalues of M.c are in 





' -\h Ell vi 
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~h v 3 


-ei" 
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Table 2 

Counting of free parameters in MSSM 



Parameters 


Conditions 


Free Parameters 


ht, hb, h T 


m w , m t 


tan (3 


v d , v u ,M 1/2 


nib, to t 


2 Extra 


mo, A A* 


ti=0, £ = 1,2 


(e.g. m h , m A ) 


Total = 9 


Total = 6 


Total = 3 


Table 3 






Counting of free parameters 


m our model 


Parameters 


Conditions 


Free Parameters 


ht, hb, h T 


mw, m t 


tan/3, ti 


Vd, V u , Mi/2 


mb, m T 




too, A, (i 


U = 


2 Extra 


Vi, ei 


(i = l,...,5) 


(e.g. m h , m A ) 


Total = 15 


Total = 9 


Total = 6 



agreement with the experimental masses of the 
leptons. After running the RGE we have a com- 
plete set of parameters, Yukawa couplings and 
soft-breaking masses mf (RGE) to study the min- 
imization. This is done by the following method: 
we solve the minimization equations for the soft 
masses squared. This is easy because those equa- 
tions are linear on the soft masses squared. The 
values obtained in this way, that we call to 2 are 
not equal to the values mf (RGE) that we got via 
RGE. To achieve equality we define a function 



rj = max 



l (RGE)\ 



m}(RGEy mj 
with the obvious property that 

V > 1 



Vi 



(41) 



(42) 



Then we adjust the parameters to minimize n. 

Before we end this section let us discuss the 
counting of free parameters in this model and 
in the minimal N=l supergravity unified version 
of the MSSM. In Table || we show this counting 
for the MSSM and in Table | for the e-model. 
Finally, we note that in either case, the sign of 
the mixing parameter /i is physical and has to be 
taken into account. 



G 



5.2. Gauge and Yukawa Unification in the 
e model 

There is a strong motivation to consider GUT 
theories where both gauge and Yukawa unification 
can achieved. This is because besides achieving 
gauge coupling unification, GUT theories also re- 
duce the number of free parameters in the Yukawa 
sector and this is normally a desirable feature. 
The situation with respect to the MSSM can be 
summarized as follows: 

• In SU(5) models, hi, = h T at Mqut- The 
predicted ratio mb/m T at Mweak agrees 
with experiments. 

• A relation between m top and tan (3 is pre- 
dicted. Two solutions are possible: low and 
high tan (3 . 

• In 50(10) and Eq models h t = ht = h T at 
Mgut- In this case, only the large tan/3 
solution survives. 

We have shown Jll| that the e-model allows b — r 
Yukawa unification for any value of tan (3 and sat- 
isfying perturbativity of the couplings. We also 
find the t — b — r Yukawa unification easier to 
achieve than in the MSSM, occurring in a wider 
high tan [3 region. 

6. Tree Level Neutrino Masses and Mix- 
ings 

6.1. Neutral fermion mass matrix 

In the basis 

V> 0T = HA', -zA 3 , Hi Hi u e , v T ) (43) 

the neutral fermions mass terms in the Lagran- 
gian are given by 

1 



C m = -^{rVM N ^ + h.c. (44) 
where the neutralino/neutrino mass matrix is 



M N 
with 



M x o ui 
m 



(45) 





\g'v u 



M 2 



-\g'vd 
\gvd 



-\gv u 



\gvd 
-hgv u 



(46) 



m 



\g'vi 



\gvi ei 



)g'v2 \gv2 o £ 2 



2i 

2 H "2 2* 

\g'n \gv?, e 3 

The mass matrix Mjv is diagonalized by 
A/"* M w AT 1 



diag(m x o, m Uj ) 



(47) 



(48) 



where (i = 1, • • • , 4) for the neutralinos, and (j = 
1, • • • , 3) for the neutrinos. 

6.2. Approximate diagonalization of mass 
matrices 

If the ffip parameters are small it is convenient 
to define |2(| the matrix 

£ = m • M-Z (49) 

If the elements of this matrix satisfy 

« 1 (50) 

then one can find an approximate solution for 
mixing matrix J\f. Explicitly we have 



& = 
62 = 

£i3 = 

£i4 = 
where 



■A, 



g'M 2 fi 
2det(M x o 

_ ghhn 
2det(M x o) 



A, 



{g 2 M 1 + g' 2 M 2 )v 2 
Adet(M x o) 



A; 



(g 2 M 1 + g^M^V! 
4det(M x o) 



A; 



A, = fivi + v d Ei 



(51) 



(52) 



From Eq. (|5l]) and Eq. ( |52| ) one can see that 
f = in the MSSM limit where e; = 0, Vi = 0. In 
leading order in £ the mixing matrix J\f is given 
by, 



A"*: 



N* 

v, T 



-e 1 - 5« + 



(53) 



The second matrix above block-diagonalizes Mn 
approximately to the form diag(m e yy, M x o) 



m e ff = - to • M x om T 



7 



M^ + M 2 g' 2 / k \ K & ^„ 

fTTFTT^ A - A m A l A M A r 54 
4det(M x o) [ AeA » ^ ^ I 



P"<7 



The sub-matrices AT and V„ in cq. diagonalize 
A4 x o and m e // 



N*M x oN^ = diag(ra x o), 

Vjm e f f V v = diag(0, 0, m„) 
where 

m„ = Tr(m e ff 



(55) 
(56) 



Ml52 + M25 ' 2 |A| 2 . (57) 



4det(7W x o) 

For V v we have ( we can rotate away one angle) 



V v = 



1 









cosf 




— sin 


#23 


sin 6 


23 


cos 6 


23 


COS 9i3 





— sin 


#13 





1 







sin 6*i3 





cos 6 


'13 



where the mixing angles can be expre: 
of the alignment vector A as follows: 



(58) 

;les can be expressed in terms 
(59) 
(60) 



tan = = 

(A*+A?)* 

tan 6*23 = 

7. One Loop Neutrino Masses and Mixings 

7.1. Definition 

The Self-Energy for the neutralino/ncutrino 



is 




j/p[p L £6.+P fl £g 
~PJ\% + P R n*] J(61) 



Then 

< c = M^( m ) + AM ;j 
where 

AM ij = [| ( n ij( m i ) + n ij( m j)) 



(62) 



p-g 



1 P J 1 P J 



p.' \ 



p-q 



1 1 

\ / 



q : \ q q \ q 
1 p j 1 J 

Figure 2. Diagrams contributing 



-\ ( m x° E ^( m i)+ m x° S ij( m j)) 
where 

s v = i(s L + s K ) n y = |(n L + n*) (64) 

and 

2 



i6 3) 



A = 



4-d 



— 7B + In 47T 



(65) 



7.2. Diagrams Contributing 

In Fig. (2) are shown the classes of diagrams 
contributing to the self-energy at one loop. These 
diagrams can be calculated in a straightforward 
way. For instance the W diagram in the £ = 1 
gauge gives 



i 

ij ~ _ 16tt2 2^, 1 \ U Ljk U Lki + U Rjk U Rki) 
k=l 

2 „,2 2 



Bi{p , m k, m w) 



n v = 

11 16tt 2 



1 5 

2^' I \ U Ljk u Rki + u Rjk u Lki) 



k=l 



m k B (p ,m k ,m w ) (66) 



where Bo and B\ are the Passarino-Veltman func- 
tions, and the coupling matrices appear in the 
vertices in the following way 
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H 



q 

-*- 



igHCJC. 




,1^ (Ol™P L + O^P R )(67) 



W, i^{OWiPL+O c ^P R ) (68) 



7.3. Gauge Invariance 

When calculating the self-energies the question 
of gauge invariance arises. In the gauge the 
sets of diagrams of Figs. |^-|| depend on £. We 
have shown that the gauge dependence cancels 
among the diagrams in each set. So in the actual 




c. 



1 Shww 



lgHC + C + 

q 



igHC.C. 



H< 



i J 

Figure 5. Set 3 



calculations we considered the tadpoles needed in 
those sets to ensure gauge invariance. The other 
tadpoles were included in the minimization pro- 
cedure of Eq. ([mJ). This is a gauge invariant 
splitting. 

7.4. The One Loop Mass Matrix 

The one-loop corrected mass matrix is given 

by 

M 1L = M° d t g + AM 1L 
where 

Mf ag =NM N N T 

Now we diagonalize the 1-loop mass matrix 
Mli g =N'M 1L N' T 

Then the mass eigenstates are related to the weak 
basis states by 



(69) 
(70) 
(71) 



mass \f^L ^,weak 

AO — - /V ?q AO 



(72) 
(73) 



with 

TV 1L = TV' TV 

The usual convention in neutrino physics 

v a = U ak v k (74) 

is recovered in our notation as 

U ak =M£ kA+a (75) 

7.5. Solar and Atmospheric Neutrino Pa- 
rameters 

Assuming hierarchy in the masses m U2 and m„ 3 
the survival probabilities for the solar and atmo- 
spheric neutrinos are 



P e = l-4C/ 2 lC / e 2 2 sin 2 )-2U^(l - Ul 



^ = 1-4^3(1-^3) sin 2 



Am 21 i 



(76) 



As [43 has to be small we neglect it and write the 
usual two neutrino mixing angle as 



sin 2 (20 12 )=4C/ e 2 1 t/ 2 2 
and 

sin 2 (2M=4[/ 2 3 (l-[/ 2 3 ) 



(77) 
(78) 
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jr. 




0.01 
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A M /A r 

Figure 6. The atmospheric angle as function of 
A M /A r for |e*| = e and A e = 0.1A T . e 2 /A has an 
upper cut of e 2 /A < 0.1 in this plot, since larger 
values lead to larger scatter for very small A p / A T 



7.6. Our Preliminary Results 

We have found {§, that if e 2 /A < 100 then the 
approximate formulas hold 



U. 



<?3 



sin tan 



sin tan 



Then if we take 
Ai < A 2 ~ A 3 





hA 2 


A?- 


hA 2 


AS- 


hA 2 



(79) 



(80) 



we immediately get maximal mixing for the at- 
mospheric neutrinos. This is shown in Figure |^ 
where we see that maximality of the mixing is 
only possible for A M = A T . 

To get bimaximality we have to fix the solar 
angle. We have discovered that if e e ~ e M ~ e T 
and A e « ~ A T we get bimaximality if the 
following sign condition applies 

(e^/er) x (A^/Ar) < (81) 

This is illustrated in Figure 0. In practice we do 
not need perfect maximality. We took 

4/5 < \A^/A T \ < 5/4 and |A e /A M | « 1 (82) 



0.001 0.01 0.1 

If 1 1 — 



•v 



0.01 0.1 



g.oo 



Figure 7. The solar angle as function of fie/e^ for 
e M = e T and A M = A r applying the sign condition. 



to fix the atmospheric angle and 

(e M /&r) x (A M /A T ) < and 0.6 < e e /e^ < 1.2(83) 

to fix the solar angle. Next we have to fix the 
masses to solve the atmospheric and solar neu- 
trino problems. We found p|j20|] that the range 

0.03GeV 2 < |A| < 0.25GeV 2 (84) 

fixes the tree level mass to reproduce the atmo- 
spheric neutrino problem. This is illustrated in 
Figu re |s| where, besides the conditions in Eqs.d83) 
and ( 83 ) to fix the angles and condition Eq. (gij) , 
all the other parameters were chosen randomly. 
Consistency of the parameters was required in the 
sense that minimization the scalar potential in- 
cluding the tadpoles was performed as well as the 
matching with the RGE solutions with universal- 
ity at GUT scale. 

Finally we have to check if it is possible to have 
masses in the range to solve the solar neutrino 
problem. We found Q that the relevant parame- 
ter for this purpose is 



(4 



e 2 )/|A|^ e 2 /|A| 



(85) 



Depending on this quantity in the range 0.01 < 
e 2 /|A| < 1 the solar neutrino problem can be 
solved: Low values give just-so solutions, high 
values tend to give large angle MSW (LA-MSW). 
This is illustrated in Figure || and in Figure |l(| 
Another question of relevance that we addressed 
was the study of the decay length of the lightest 
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). All points 
obey Eqs. (|2|), (|8j and (|8j) no further cut ap- 
plied except that 0.3 < e 2 /|A| < 1. 



Figure 8. Am^ versus sin 2 (2^ 
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Figure 9. Ato 2 2 as a function of e 2 /|A|. A box is 
drawn to guide the eye. 
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neutralino. This is important, because if the de- 
cay length is greater than the detector, then in 
practice it will be invisible like in the MSSM. As 
we can see in Figure [lIJ that is not the case, the 
neutralino decays well inside the detector leading 
to novel signatures. 
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Figure 10. Ato 2 2 versus sin 2 (29 so i) for those 
points which have Am 2 3 and sin 2 (2# atm ) correct 
and at the same time fit either the vacuum or 
LA-MSW solutions. 



There is a viable model for SBRP that leads to 
a very rich phenomenology, both at laboratory 
experiments, and at present (LEP) and future 
(LHC, LNC) accelerators. We have shown that 
the radiative breaking of both the Gauge Sym- 
metry and Rp can be achieved. In these type of 
models neutrinos have mass and can decay thus 
avoiding the critical density argument. They also 
can evade the BBN limits on a v r on the MeV 
scale. Most of these phenomenology can be de- 
scribed by an effective model with bilinear ex- 
plicit Rp violation. We have calculated the one- 
loop corrected masses and mixings for the neu- 
trinos in a completely consistent way, including 
the RG equations and correctly minimizing the 
potential. We have shown that it is possible to 
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get bimaximal solutions for both the atmospheric 
and solar neutrino problems. We emphasize that 
the lightest neutralino decays inside the detectors, 
thus leading to a very different phenomenology 
than the MSSM. 
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